function [mu,w]=S_1D(x)
% given order x, returns the directions and weights

mu_table(1,2) = -0.5773502691;
mu_table(2,2) = -mu_table(1,2);

w_table(1,2)  = 1.0;
w_table(2,2)  = w_table(1,2);

mu_table(1,4) = -0.8611363115;
mu_table(2,4) = -0.3399810435;
mu_table(3,4) = -mu_table(2,4);
mu_table(4,4) = -mu_table(1,4);

w_table(1,4)  = 0.3478548451;
w_table(2,4)  = 0.6521451549;
w_table(3,4)  = w_table(2,4);
w_table(4,4)  = w_table(1,4);

mu_table(1,8) = -0.9602898564;
mu_table(2,8) = -0.7966664774;
mu_table(3,8) = -0.5255324099;
mu_table(4,8) = -0.1834346424;
mu_table(5,8) = -mu_table(4,8);
mu_table(6,8) = -mu_table(3,8);
mu_table(7,8) = -mu_table(2,8);
mu_table(8,8) = -mu_table(1,8);

w_table(1,8)  = 0.1012285363;
w_table(2,8)  = 0.2223810344;
w_table(3,8)  = 0.3137066459;
w_table(4,8)  = 0.3626837834;
w_table(5,8)  = w_table(4,8);
w_table(6,8)  = w_table(3,8);
w_table(7,8)  = w_table(2,8);
w_table(8,8)  = w_table(1,8);

mu_table(1,16)  = -0.9894009349916499;
mu_table(2,16)  = -0.9445750230732326;
mu_table(3,16)  = -0.8656312023878318;
mu_table(4,16)  = -0.755404408355003;
mu_table(5,16)  = -0.6178762444026438;
mu_table(6,16)  = -0.4580167776572274;
mu_table(7,16)  = -0.2816035507792589;
mu_table(8,16)  = -0.09501250983763744;
mu_table(9,16)  = -mu_table(8,16);
mu_table(10,16) = -mu_table(7,16);
mu_table(11,16) = -mu_table(6,16);
mu_table(12,16) = -mu_table(5,16);
mu_table(13,16) = -mu_table(4,16);
mu_table(14,16) = -mu_table(3,16);
mu_table(15,16) = -mu_table(2,16);
mu_table(16,16) = -mu_table(1,16);

w_table(1,16)  = 0.027152459411742137;
w_table(2,16)  = 0.06225352393877526;
w_table(3,16)  = 0.09515851168244673;
w_table(4,16)  = 0.12462897125537288;
w_table(5,16)  = 0.1495959888165501;
w_table(6,16)  = 0.16915651939499357;
w_table(7,16)  = 0.1826034150449233;
w_table(8,16)  = 0.1894506104550685;
w_table(9,16)  = w_table(8,16);
w_table(10,16) = w_table(7,16);
w_table(11,16) = w_table(6,16);
w_table(12,16) = w_table(5,16);
w_table(13,16) = w_table(4,16);
w_table(14,16) = w_table(3,16);
w_table(15,16) = w_table(2,16);
w_table(16,16) = w_table(1,16);

mu_table(1,32)  = -0.9972638618494816;
mu_table(2,32)  = -0.9856115115452684;
mu_table(3,32)  = -0.9647622555875064;
mu_table(4,32)  = -0.9349060759377397;
mu_table(5,32)  = -0.8963211557660521;
mu_table(6,32)  = -0.84936761373257;
mu_table(7,32)  = -0.7944837959679424;
mu_table(8,32)  = -0.7321821187402897;
mu_table(9,32)  = -0.6630442669302152;
mu_table(10,32) = -0.5877157572407623;
mu_table(11,32) = -0.5068999089322294;
mu_table(12,32) = -0.42135127613063533;
mu_table(13,32) = -0.33186860228212767;
mu_table(14,32) = -0.2392873622521371;
mu_table(15,32) = -0.1444719615827965;
mu_table(16,32) = -0.04830766568773832;
mu_table(17,32) = -mu_table(16,32);
mu_table(18,32) = -mu_table(15,32);
mu_table(19,32) = -mu_table(14,32);
mu_table(20,32) = -mu_table(13,32);
mu_table(21,32) = -mu_table(12,32);
mu_table(22,32) = -mu_table(11,32);
mu_table(23,32) = -mu_table(10,32);
mu_table(24,32) = -mu_table(9,32);
mu_table(25,32) = -mu_table(8,32);
mu_table(26,32) = -mu_table(7,32);
mu_table(27,32) = -mu_table(6,32);
mu_table(28,32) = -mu_table(5,32);
mu_table(29,32) = -mu_table(4,32);
mu_table(30,32) = -mu_table(3,32);
mu_table(31,32) = -mu_table(2,32);
mu_table(32,32) = -mu_table(1,32);

w_table(1,32)  = 0.007018601052751814;
w_table(2,32)  = 0.016274489134646575;
w_table(3,32)  = 0.025392016655236452;
w_table(4,32)  = 0.03427396600007983;
w_table(5,32)  = 0.04283582157141583;
w_table(6,32)  = 0.05099805598882692;
w_table(7,32)  = 0.058684062962362976;
w_table(8,32)  = 0.06582222058112427;
w_table(9,32)  = 0.07234579320304929;
w_table(10,32) = 0.07819389600818424;
w_table(11,32) = 0.08331192422405545;
w_table(12,32) = 0.08765209300583078;
w_table(13,32) = 0.09117387869558631;
w_table(14,32) = 0.093844399080795;
w_table(15,32) = 0.095638720079275;
w_table(16,32) = 0.0965400885147278;
w_table(17,32) = w_table(16,32);
w_table(18,32) = w_table(15,32);
w_table(19,32) = w_table(14,32);
w_table(20,32) = w_table(13,32);
w_table(21,32) = w_table(12,32);
w_table(22,32) = w_table(11,32);
w_table(23,32) = w_table(10,32);
w_table(24,32) = w_table(9,32);
w_table(25,32) = w_table(8,32);
w_table(26,32) = w_table(7,32);
w_table(27,32) = w_table(6,32);
w_table(28,32) = w_table(5,32);
w_table(29,32) = w_table(4,32);
w_table(30,32) = w_table(3,32);
w_table(31,32) = w_table(2,32);
w_table(32,32) = w_table(1,32);

mu = mu_table(1:x,x);
w  = w_table(1:x,x);
